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Abstract.  In  the  DLR-  high  vacuum  plume  test  facility  STG  experiments  have  been  performed  with  the 
aim  to  validate  the  available  plume  flow  model.  For  this  purpose  fully  expanding  simulated  plumes  with  N2 
as  test  gas  have  been  chosen  as  simplest  ceise  (no  chemical  reactions,  constant  ratio  of  specific  heats).  The 
essential  parameter  to  be  varied  was  the  nozzle  flow  Reynolds  number,  the  quantities  to  be  measured  were 
the  Pitot  pressure  at  the  nozzle  exit  and  the  molecule  number  flux  in  the  plume.  It  turned  out  that  a  new 
model  had  to  be  written  taking  into  account  in  a  more  appropriate  way  the  development  of  the  boundary 
layer  in  the  nozzle  and  the  initial  plume  spreading.  The  result  is  a  code  which  describes  the  plume  flow  in 
satisfying  agreement  with  the  experimental  data. 


INTRODUCTION 

In  the  past  many  attempts  have  been  made  to  describe  the  fully  expanding  thruster  plume  flow  by  relatively 
simple  models  as  well  as  by  sophisticated  calculations,  but  only  a  few  experimental  investigations  have  been 
performed  in  the  rarefled  regime,  which  is  by  far  the  larger  and  more  important  part  of  this  flow  fleld  with 
respect  to  impingement  problems  [1] .  At  the  last  RGD  we  have  reported  on  molecule  number  flux  measurements 
in  fully  expanding  free  jets  and  plumes  with  test  gases  N2  and  H2  in  the  DLR-high  vacuum  plume  test  facility 
STG  [2].  These  experimental  investigations  have  been  continued  systematically  with  the  aim  to  validate  the 
available  DLR-model  [3]  or  to  improve  it.  Since  there  are  various  parameters  influencing  an  actual  thruster 
plume  flow  (nozzle  shape,  chemical  reactions  and  condensation,  variation  of  the  ratio  of  specific  heats  with 
decreasing  temperature  during  the  expansion,  et  cetera)  it  is  appropriate  to  reduce  their  number  at  first  and  to 
choose  the  experimental  conditions  as  simple  as  possible.  In  this  paper  we  will  first  present  a  comparison  of  the 
results  of  such  experiments  and  of  calculations  with  the  DLR-model  [3].  Since  the  agreement  turned  out  not 
to  be  satisfactory,  an  analysis  of  the  reasons  follows.  Based  on  this  a  new  model  is  developed  and  introduced. 


EXPERIMENTAL  RESULTS 

The  test  facility  STG  and  the  free  molecule  pressure  probe  to  measure  the  molecule  number  flux  have  been 
introduced  previously  [4].  In  our  tests  the  nozzle  shape  was  conical  (DAS A  0.5N  monopropellant  hydrazine 
thruster  nozzle:  Half  angle  15°,  length  Ie  =  7.7bmm,  exit  diameter  ds  =  4.75mm,  throat  diameter  d*  = 
0.6  mm)  and  Nitrogen  was  selected  as  test  gas  to  have  a  single  component  gas  flow  without  interactions  and 
phenomena  possibly  occuring  in  a  gas  mixture.  The  stagnation  temperatures  To  were  300  K  and  600  if;  thus 
a  constant  ratio  of  specific  heats  k  was  maintained  throughout  the  flow.  Then  the  stagnation  pressure  was 
varied  providing  a  variation  of  the  nozzle  boundary  layer  thickness  5  ~  =ue  Pe  IeIp{To):  nozzle 

flow  Reynolds  number;  pe,  ue  are  the  velocity  and  density  at  the  nozzle  exit,  calculated  with  the  assumption 
of  a  pure  isentropic  flow,  and  p  is  the  viscosity).  Using  a  polar  coordinate  system  r,  6  (origin  at  the  nozzle  exit 
center,  9  =  0°  corresponds  to  the  plume  axis),  an  angular  number  flux  profile  family,  measured  at  fixed  radius 
r  =  0.5  m  in  the  range  —90°  <6<  90°  with  a  variation  of  Rce,  is  shown  in  Fig.  1  together  with  the  calculated 
results  for  Rce  =  800  and  3200.  Different  attempts  were  made  to  get  a  better  agreement,  for  example  by 
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FIGURE  1.  Angular  profiles  of  the  number  flux  hi  with  diflferent  Res',  calculations  with  the  former  DLR-model 

using  other  simple  functions  (up  to  now  a  cos-function  was  used  for  the  isentropic  core  and  an  exponential 
decay  to  describe  the  number  flux  distribution  in  the  boundary  layer  expansion  region).  But  it  turned  out  that 
two  main  features  of  the  profile  family  in  Fig.  1  could  not  be  reproduced  in  this  way,  namely  the  inflection 
points  in  the  vicinity  of  6  =  ±40°  and  the  relatively  increasing  concentration  of  the  number  flux  in  the  near 
axis  region  with  increasing  nozzle  Reynolds  number.  Hence  it  became  clear  that  the  proper  description  of  the 
nozzle  boundary  layer  and  of  its  influence  on  the  isentropic  core  is  decisive  and  it  turned  out  that  the  concept 
used  up  to  now  is  too  simple. 


MODELLING 

The  isentropic  core  flow  in  the  nozzle  was  treated  as  a  one-dimensional  diverging  streamtube  which  is  sur¬ 
rounded  by  a  boundary  layer  (see  Fig.  2).  The  calculation  procedure  starts  in  the  throat  with  the  assumption 
of  a  pure  isentropic  flow  along  the  length  Ax  in  flow  direction.  For  the  next  interval  Ax  downstream  the 
isentropic  flow  conditions  have  then  to  be  matched  with  the  boundary  layer  conditions  at  the  wall  in  this 
interval  and  with  the  following  functions  describing  the  velocity  and  temperature  in  the  boundary  layer,  which 
is  assumed  to  be  hypersonic  and  laminar:  Hantzsche  and  Wendt  [5,6]  showed  in  their  calculations,  that  the 
velocity  profile  u{y)  of  the  plane  hypersonic  laminar  boundary  layer  could  be  approximated  satisfactorily  by  a 
linear  function 

y 

u{y)=u^-  with  0<y<6, 

0 

where  y  denotes  the  distance  from  the  wall,  Uoo  the  velocity  of  the  flow  outside  the  boundary  layer  and  S 
the  boundary  layer  thickness.  In  [7]  Hanztsche  and  Wendt  demonstrate  that  the  boundary  layer  equations 
for  a  flat  plate  could  be  transformed  into  the  coordinates  of  a  cone  by  a  simple  constant  (C  =  0.577).  They 
demonstrate  this  for  an  exterior  flow,  and  we  adopt  this  transformation  for  the  (interior)  nozzle  flow.  The  wall 
temperature  T^,  which  has  been  determined  by  Hantzsche  and  Wendt  [7],  Korkegi  [8]  and  van  Driest  [9]  with 
the  assumption  of  an  adiabatic  wall,  could  be  approximated  by  the  recovery  temperature: 

Tw  =  Too  ^1  ± 

So  the  temperature  profile  T{y)  could  be  written  as 

T{y)  =  (Tu)  -  Too)  cos^'®  with  0<y  <St, 

where  St  denotes  the  temperature  boundary  layer  thickness  and  Too  the  free  stream  temperature.  The  exponent 
of  the  cos-function  has  been  determined  empirically  by  fitting  based  on  Pitot  pressure  measurements  (see  Fig. 
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4).  In  this  calculation  procedure  the  nozzle  flow  Reynolds  number  is  Reg  =  (s  denotes  the  distance 

along  the  nozzle  wall  and  J  =  C  )•  With  the  Prandtl  number  Pr  =  0.74  for  N2  the  temperature  boundary 

layer  thickness  5t  could  be  estimated  with  ^  fs  '/Pr.  So  the  displacement  thickness  (5i  of  this  intervall  Ax, 
which  reduces  the  cross-section  of  the  isentropic  core  in  the  nozzle,  can  be  calculated  with  the  integration  up 
to  5t-  Matching  of  isentropic  core  and  boundary  layer  in  every  interval  is  achieved  by  iteration.  The  values  of 
the  former  interval  are  used  as  initial  values  for  the  iteration  of  the  next  interval.  So  the  calculation  marches 
downstream  to  the  nozzle  exit  (see  Fig.  2).  This  concept  has  also  been  used  by  Lengrand  [10],  but  he  integrates 
the  von  Kdrmdn  integral  equation  to  determine  the  displacement  thickness. 


nozzle  throat 


nozzle  exit 


FIGURE  2.  Sketch  to  explain  the  downstream  marching  iteration  of  the  nozzle  flow 


Having  calculated  the  complete  nozzle  flow  in  this  way  (example  see  Fig.  3),  Pitot  pressure  profiles  Pt2{y)  at 
the  nozzle  exit  were  determined  both  by  calculation  and  measurement  and  sketched  in  Fig.  4  (the  details  of  the 
measurement,  of  the  concept  and  of  the  calculation  of  the  nozzle  flow  are  outlined  in  [11]).  The  agreement  is 
considered  as  satisfactory  and  the  result  of  this  simple  approximation  of  the  nozzle  flow  was  therefore  considered 
as  suitable  to  continue  with  the  plume  expansion  modelling. 

In  our  previous  model  [3]  straight  streamlines  originated  directly  at  the  nozzle  exit  with  an  asssumed  proper 
assignment  of  their  characterising  angle  6  with  the  plume  axis  to  the  exit  coordinate  y  perpendicular  to  the 
plume  axis.  Knowing  now  the  number  flux  distribution  at  the  nozzle  exit  (from  the  calculation)  and  in  the 
plume  (from  the  experiments),  it  is  of  course  possible  to  find  a  better,  fitting  assignment.  However,  the  former 
concept  means  physically  that  the  expansion  potential  of  the  gas  downstream  of  the  nozzle  is  neglected  (this 
is  clearly  to  be  seen  in  Fig.  1).  To  consider  this  expansion  potential  we  tried  to  apply  a  concept  to  the  flow 
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FIGURE  4.  Normalized  Pitot  pressure  profiles  at  the  nozzle  exit;  po:  Stagnation  pressure 


FIGURE  5.  Sketch  of  the  surface  elements  at  the  nozzle  exit  AAe  and  on  the  reference  half  sphere  A^s 


conditions  at  the  nozzle  exit,  which  is  as  drastic  as  the  one  of  the  freezing  surface  used  to  describe  the  transition 
to  free  molecule  flow,  i.e.  a  superposition  of  contributions  from  different  annular  surface  elements  at  the  nozzle 
exit  (Fig.  5).  The  justification  at  first  is  that  the  flow  near  the  nozzle  lip,  at  least  immediately  downstream  in 
the  off-axis  region  is  rarefied,  and  the  thermal  motion  of  the  gas  provides  a  lateral  spreading  of  the  molecules 
with  crossing  particle  paths. 

The  distribution  of  the  gas  portions  originating  from  a  surface  element  AA^  (Fig.  5)  into  the  different 
directions  ip  is  described  by  an  empirical  weight  function.  This  weight  function  is  a  Gaufi  distribution  which 
is  defined  as 

w{ip)  =  . 

For  the  determination  of  a  and  b  we  first  use  the  mass  conservation  which  is  fulfilled  by  itself  (see  Fig.  6): 

2tv  7r/2 

Ip  w (ip)  dip  dip  =  1  and  this  leads  to  a  = 

0  0 
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FIGURE  6.  Weight  function  w(^)  to  demonstrate  the  geometric  relations 


The  quantity  b  is  set  to  6  =  It  depends  on  the  ratio  of  specific  heats  k  and  the  Mach  number 

Ma,  since  they  are  the  most  important  parameters  which  determine  the  expansion  potential  at  a  certain  surface 
element  AAe  at  the  nozzle  exit.  The  resulting  function  w(^)  for  k  =  1.4  with  Ma  as  parameter  is  shown  in 
Fig.  7.  This  particular  introduction  of  b  was  empirically:  In  the  isentropic  core  it  was  found  that  with  the 
calculated  Mach  numbers  at  the  nozzle  exit  {MaE,is  ^  5.5)  most  of  the  mass  flux  is  spread  into  a  region  of 
0  <  0  <  30°  (see  also  Fig.  1),  and  if  we  consider  the  continuum  near  axis  flow  at  the  nozzle  exit  as  single 
stream  tube,  the  function  applied  to  this  flow  yields  essentially  hi{6)  as  it  is  found  in  Fig.  1  in  the  plume 
near  axis  region.  The  application  to  surface  elements  AAe  delivers  a  refinement  but  no  basic  change.  This  is 
the  reason  why  we  can  apply  this  concept  with  the  weight  function  to  the  whole  nozzle  exit  flow. 

To  simplify  the  calculation  the  difference  between  6  and  ^p  (Fig.  5)  for  the  same  surface  element  AAg  is 
neglected:  6  =  ip,  and  the  number  of  the  surface  elements  at  the  nozzle  exit  and  at  the  half  sphere  r^e/  is  set 
to  k,  which  leads  to 

T*  E 

Aro  =  -y  ,  rD,i  =  i  Am  ,  for  i  =  0, . . . ,  fc  -  1 

k 


iP 

FIGURE  7.  Weight  function  w{ip)  for  different  Mach  numbers  and  k  =  1.4 
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Hence  the  number  flux  ni{rref,Oj)  through  a  surface  element 

2  ) 

on  a  half  sphere  at  the  distance  Vref  can  be  calculated  as  the  sum  of  all  weigthed  number  fluxes  from  the 
surface  elements 


COS 


8j 


AAE,i  =  n  (Ar|,  +  2rD,i  Arn) 


at  the  nozzle  exit  : 


1 

AAgj 


f  \  /  \  AAE,i 
2^  WE,ij  u{yi)  p(yi)—— 

i=0 


,  with 


Vi  =  rE- 


( 


rD,i  + 


ArD\ 

2  ) 


The  values  of  the  weight  function  are  given  by 


WE,ij 


,  with  Oj,  hi  parameters  of  w{'ip)  depending  on  Ma{yi). 


To  avoid  freezing  on  the  outer  surface  elements  AAg  during  the  superposition  from  the  nozzle  exit  plane  to 
the  half  sphere  r^e/ ,  the  calculation  of  the  number  flux  is  repeated  with  decreasing  radius  r^e/  until  the  flow 
on  no  outer  surface  element  is  frozen.  Starting  at  this  first  half  sphere  r^ef ,  the  calculation  of  the  whole  plume 
marches  downstream  with  an  assumed  density  decay  p  ~  like  in  the  former  model  (example  see  Fig. 
8;  x,y  denote  in  this  and  the  next  figure  cartesian  coordinates  with  the  origine  in  the  nozzle  exit).  In  this 
region  freezing  is  taken  into  consideration  by  the  calculation  of  a  freezing  surface  [3]  at  P/  =  2  {Bird’s  freezing 
parameter;  see  Fig.  9). 

Two  results  of  the  present  calculation  in  comparison  to  the  experiments  are  shown  in  Fig.  10  with  the 
experimental  curves  taken  from  Fig.  1.  The  Reynolds  numbers  represent  two  extreme  cases  of  boundary  layer 


FIGURE  9.  Location  of  the  calculated  freezing  surface  with  different  Res 
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FIGURE  10.  Results  of  the  present  calculation  in  comparison  to  the  experiments 


FIGURE  11. 

DLR-model 


thickness  regarding  the  actual  (monopropellant  hydrcizine)  thruster  operation.  Qualitatively  the  curves  have 
the  same  features,  especially  the  inflection  points  at  0  w  40“  are  reproduced.  The  agreement  on  the  axis  is 
very  good,  in  the  off-axis  region  the  discrepancy  has  been  reduced  drastically  compared  to  that  of  our  previous 
model  (see  Fig.  1). 

Just  for  completeness  we  want  to  mention  that  in  addition  to  the  number  flux  other  flow  quantities  like  the 
velocity  u  and  the  temperature  T  can  be  calculated.  During  the  superposition  of  the  weighted  contributions  the 
other  flow  quantities  are  weighted  by  their  kinetic  energy,  and  they  are  calculated  by  the  method  of  effective 
stagnation  conditions  (see  [3]).  A  direct  comparison  with  experimental  results  is  not  yet  possible,  but  from  the 
physical  plausibility  it  seems  that  the  results  obtained  with  the  new  model  are  more  realistic  (see  Fig.  11). 
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CONCLUSIONS 


The  aim  of  our  study  was  to  validate  or  to  improve  the  DLR-  plume  model  based  on  experiments  with 
simulated  plumes  in  STG.  The  measurements  especially  in  the  boundary  layer  expansion  region  showed  that 
the  available  model  concept  was  not  sufficient,  but  that  the  nozzle  flow  boundary  layer  characteristics  and  the 
initial  turning  of  the  streamlines  during  the  free  expansion  had  to  be  taken  stronger  into  the  considerations. 
We  used  available  boundary  layer  calculation  concepts  and  applied  them  to  the  nozzle  flow.  By  means  of 
Pitot  pressure  measurements  at  the  nozzle  exit  location  the  reliability  of  the  nozzle  flow  calculations  could  be 
shown.  A  weight  function  which  is  essentially  determined  by  Ma  and  k  was  then  introduced  to  perform  the 
transition  from  the  flow  at  the  nozzle  exit  to  a  reference  half  sphere  downstream.  Thus  the  initial  turning 
of  the  streamlines  in  the  free  expansion  was  modelled.  The  flow  quantities  at  this  reference  surface  are  the 
starting  condition  for  the  calculation  of  the  whole  plume.  We  found  that  the  overall  resulting  number  flux 
distribution  in  the  plume  shows  a  satisfying  agreement  with  the  measured  one,  especially  when  we  regard  that 
the  magnitude  of  this  flow  quantity  changes  within  this  local  distribution  considerably  (relative  change  from 
about  10^  at  0  =  0°  to  1  at  0  =  90°) .  A  simple  plume  description  was  maintained  (the  duration  of  the  present 
calculation  on  a  normal  personal  computer  is  about  one  minute  to  estimate  the  nozzle  flow  and  the  rarefied 
regime  up  to  a  distance  of  r  =  5  m),  and  it  is  intended  to  further  develop  the  concept  to  more  complicated 
cases,  at  first  using  H2,  where  k  is  not  constant,  and  test  gas  mixtures  with  N2  and  H2  to  simulate  more 
realistic  actual  thruster  plume  conditions. 
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